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Abstract 



Lithium niobate photonic circuits have the salutary property of permitting the generation, transmis- 
sion, and processing of photons to be accommodated on a single chip. Compact photonic circuits such 
as these, with multiple components integrated on a single chip, are crucial for efficiently implementing 
quantum information processing schemes. We present a set of basic transformations that are useful for 
manipulating modal qubits in Ti:LiNb03 photonic quantum circuits. These include the mode analyzer, 
vq ■ a device that separates the even and odd components of a state into two separate spatial paths; the 

iys | mode rotator, which rotates the state by an angle in mode space; and modal Pauli spin operators that 

effect related operations. We also describe the design of a deterministic, two-qubit, single-photon, CNOT 
gate, a key element in certain sets of universal quantum logic gates. It is implemented as a Ti:LiNb03 
photonic quantum circuit in which the polarization and mode number of a single photon serve as the 
control and target qubits, respectively. It is shown that the effects of dispersion in the CNOT circuit 
can be mitigated by augmenting it with an additional path. The performance of all of these components 
are confirmed by numerical simulations. The implementation of these transformations relies on selective 
and controllable power coupling among single- and two-mode waveguides, as well as the polarization 
sensitivity of the Pockels coefficients in LiNb03. 

X" 

1 Introduction 

We recently investigated the possibility of using spontaneous parametric down-conversion (SPDC) in two- 
mode waveguides to generate guided-wave photon pairs entangled in mode number, using a cw pump source. 
If one photon is generated in the fundamental (even) mode, the other will be in the first-order (odd) mode, and 
vice versa pQ. We also considered a number of detailed photonic-circuit designs that make use of Ti:LiNb03 
diffused channel, two-mode waveguides for generating and separating photons with various combinations 
of modal, spectral, and polarization entanglement [2j. Selective mode coupling between combinations of 
adjacent single-mode and two-mode waveguides is a key feature of these circuits. 

Although potassium titanyl phosphate (KTiOPOzi, KTP) single- and multi-mode waveguide structures 
have also been used for producing spontaneous parametric down-conversion [3HZ], it appears that only the 
generation process, which makes use of a pulsed pump source, has been incorporated on-chip. Substantial 
advances have also recently been made in the development of single-mode silica-on-silicon waveguide quantum 



circuits 015], with an eye toward quantum information processing applications [TMT5] . For these materials, 
however, the photon-generation process necessarily lies off-chip. 

Lithium niobate photonic circuits have the distinct advantage that they permit the generation, transmis- 
sion, and processing of photons all to be achieved on a single chip 0. Moreover, lithium niobate offers a 
number of ancillary advantages: 1) its properties are well-understood since it is the basis of integrated-optics 
technology |16j ; 2) circuit elements, such as two-mode waveguides and polarization-sensitive mode-separation 
structures, have low loss [2]; 3) it exhibits an electro-optic effect that can modify the refractive index at rates 
up to tens of GHz and is polarization-sensitive (TTJ Sec. 20. ID]; and 4) periodic poling of the second- 
order nonlinear optical coefficient is straightforward so that phase- matched parametric interactions [THJQIJ] , 
such as SPDC and the generation of entangled-photon pairs 0I][2T], can be readily achieved. Moreover, 
consistency between simulation and experimental measurement has been demonstrated in a whole host of 
configurations [22H26J . To enhance tolerance to fabrication errors, photonic circuits can be equipped with 
electro-optic adjustments. For example, an electro-optically switched coupler with stepped phase-mismatch 
reversal serves to maximize coupling between fabricated waveguides |27J[55] . 

Compact photonic circuits with multiple components integrated on a single chip, such as the ones con- 
sidered here, are likely to be highly important for the efficient implementation of devices in the domain of 
quantum information science. The Controlled-NOT (CNOT) gate is one such device. It plays an important 
role in quantum information processing, in no small part because it is a key element in certain sets of uni- 
versal quantum logic gates (such as CNOT plus rotation) that enable all operations possible on a quantum 
computer to be executed [TTJ[T5][25]|3n] • Two qubits are involved in its operation: a control and a target. 
The CNOT gate functions by flipping the target qubit if and only if the control qubit is in a particular state 
of the computational basis. Two separate photons, or, alternatively, two different degrees-of-freedom of the 
same photon, may be used for these two qubits. A deterministic, two-qubit, single-photon, CNOT gate 
was demonstrated using bulk optics in 2004 [31] . More recently, a probabilistic, two-photon, version of the 
CNOT gate was implemented as a silica-on-silicon photonic quantum circuit; an external bulk-optics source 
of polarization qubits was required, however 0. It is worthy of mention that qubit decoherence is likely to 
be minimal in photonic quantum circuits; however, decoherence resulting from loss in long waveguides can 
be mitigated by the use of either a qubit amplifier [32 or teleportation and error-correcting techniques 33 ■ 

This paper describes a set of basic building blocks useful for manipulating modal qubits in Ti:LiNb03 
photonic quantum circuits. Section 2 provides a brief description of the geometry and properties of the 
diffused channel Ti:LiNb03 waveguides used in the simulations. Modal qubits are characterized in Sec. 3. 
Section 4 addresses the coupling of modes between two adjacent waveguides; several special cases are high- 
lighted. The principle of operation of the mode analyzer, which separates the even and odd components of an 
incoming state into two separate spatial paths, is set forth in Sec. 5, as are the effects of the modal Pauli spin 
operator a z . The mode rotator, which rotates the state by an angle in mode space, is examined in Sec. 6, 
as is the modal Pauli spin operator a x . Section 7 is devoted to describing the design of a deterministic, 
two-qubit, single-photon, CNOT gate implemented as a Ti:LiNb03 photonic quantum circuit, in which the 
polarization and mode number of a single photon serve as the control and target qubits, respectively. The 
conclusion is presented in Sec. 8. 

2 Diffused channel Ti:LiNb03 waveguides 

All of the simulations presented in this paper refer to structures that make use of Ti:LiNb03 diffused channel 
waveguides, as illustrated in Fig. [TJ These waveguides are fabricated by diffusing a thin film of titanium 
(Ti), with thickness 5 sa 100 nm and width w, into a z-cut, y-propagating LiNb03 crystal. The diffusion 
length D is taken to be the same in the two transverse directions: D = 3 /im. The TE mode polarized in the 
cc-direction sees the ordinary refractive index n , whereas the TM mode polarized in the z-direction (along 
the optic axis) sees the extraordinary refractive index n e . 

The ordinary and extraordinary refractive indices may be calculated by making use of the Sellmeier equa- 
tions [TTJ Chap. 5] , [34|f35] . The refractive-index increase introduced by titanium indiffusion is characterized 
by An = 2Sp erf(w/2D) j ypK D, where p = 0.47 and 0.625 for n and n e , respectively [36 . To accommodate 
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Figure 1: Cross-sectional view of the fabrication of a diffused channel Ti:LiNb03 waveguide (not to scale). 
A thin film of titanium of thickness S s=a 100 nm and width w is diffused into a z-cut, y-propagating LiNb03 
crystal. The diffusion length D = 3 /im. 

wavelength dispersion, An can be modified by incorporating the weak factor £ = 0.052 + 0.065/A 2 , where 
the wavelength A is specified in /im [37] . We calculate the effective refractive index n c g of a confined mode in 
two ways: 1) by using the effective- index method described in 38 ; and 2) by making use of the commercial 
photonic and network design software package RSoft. The propagation constant of a guided mode is related 
to n e ff via f3 = 2im e s/\. 

Applying a steady electric field to this structure in the z-direction (along the optic axis) changes the 
ordinary and extraordinary refractive indices of this uniaxial (trigonal 3m) material by — -^v? r\^V j d and 
— \r? e r- i %V I d , respectively [HI Example 20.2-1], where V is the applied voltage; d is the separation between 
the electrodes; and r%3 and r^ are the tensor elements of the Pockels coefficient, which have values 10.9 and 
32.6 pm/V, respectively [35] . 

3 Modal qubits 

A qubit is a pure quantum state that resides in a two-dimensional Hilbert space. It represents a coherent 
superposition of the basis states, generally denoted |0) and |1). A qubit can be encoded in any of several 
degrees-of- freedom of a single photon, such as polarization [3S], spatial parity [20], or the mode number 
of a single photon confined to a two-mode waveguide [T][2]- The Poincare sphere provides a geometrical 
representation for the state of a modal qubit, much as it does for polarization [TT1 Sec. 6.1 A] and spatial 
parity [41"] . 

Indeed, polarization offers an intrinsically binary basis and is often used to realize a qubit. However, 
the spatial modes of a photon in a two-mode waveguide, one of which is even and the other odd, are also 
binary and can therefore also be used to represent a qubit. Modal qubits are particularly suited to photonic 
quantum circuits since they can be both generated and easily transformed on-chip by making use of elements 
such as mode analyzers, mode rotators, and two-mode electro-optic directional couplers. The modal space 
of a two-mode waveguide therefore offers an appealing alternative to polarization for representing qubits in 
quantum photonic circuits. 

The comparison between modal and spatial-parity qubits is instructive. Spatial-parity qubits are defined 
on a 2D Hilbert space in which the ID transverse spatial modes of the photon are decomposed into even 
and odd spatial-parity components [40H42] . Modal qubits also relate to parity, but in a simpler way. They 
are defined on a 2D Hilbert space in which the bases are a single ID even-parity function and a single ID 
odd-parity function. These two functions are the fundamental (even, m = 0) and first-order order (odd, 
m = 1) transverse spatial eigenmodes of the Helmholtz equation for a two- mode waveguide. 

Photon pairs can be exploited for use in quantum photonic circuits [UlEjj as well as for producing 
heralded single-photon pure states [43] in well-defined spatiotemporal modes, which are required for many 
quantum information technology applications such as quantum cryptography [44] and linear optical quantum 
computing [30]. Care must be taken, however, to ensure that the intrinsic quantum correlations between the 
twin photons are eliminated so that the surviving photon is in a pure state [45H47] . One way of achieving this 
is to generate the twin photons with a factorable joint amplitude [4"8H5~T] . We have previously shown that a 



Type-0 interaction could be used to generate photon pairs that are degenerate in frequency and polarization, 
but with opposite mode number [2j Sec. 3]. Coupling these photons into two single-mode waveguides would 
allow one of these photons to be used to herald the arrival of the other. The heralded photon could then 
be coupled into a two-mode waveguide which, with the addition of a mode rotator, would serve as a source 
of modal qubits. Such a source would be analogous to the one fashioned from bulk optics by Fiorentino et 
al. [3T] using Type-II SPDC. However, the Type-0 source of modal qubits described above would be on-chip 
and would also make use of the strongest nonlinear component of the second-order tensor, CI33, thereby 
enhancing the efficiency of the interaction [52] . 

The quantum state of a single photon in a two-mode waveguide, assuming that its polarization is TE or 
TM, can be expressed as |\&) = a\\e) + CK2I0), where |e) and \o) represent the even and odd basis states, 
respectively; and a\ and a 2 are their weights. All operations on the single-photon state are effected via 
auxiliary adjacent waveguides, which are sometimes single-mode and sometimes two-mode. We exploit the 
concepts of selective and controllable coupling between waveguides, together with the isomorphism between 
waveguide coupling and the SO(2) rotation matrix, to design a mode analyzer, a mode rotator, modal Pauli 
spin operators, and a CNOT gate useful for quantum information processing. 



4 Mode coupling between adjacent waveguides 

The coupling between two lossless, single-mode waveguides is described by a unitary matrix T that takes 
the form [HI Sec. 8.5B] 

A -jB 
-jB* A* 



(1) 



where A = exp(j A/3 L/2) [cos jL - j (A/3/27) sin 7-£] and B = (k/j) exp(jA/3L/2) sinjL . Here, A/3 is 
the phase mismatch per unit length between the two coupled modes; L is the coupling interaction length; n 
is the coupling coefficient, which depends on the widths of the waveguides and their separation as well as on 



the mode profiles; 7 



7 A/3 2 ; and the symbol * represents complex conjugation. 



This unitary matrix T can equivalently be written in polar notation as 



cos (6/2) exp(j(f> A ) 
-j sin (0/2) exp(-j</>s) 



-j sin (6/2) exp(.7</!> B ) 
cos (9/2) exp (-j4>a) 



(2) 



where 9 = 2sin _ [(k/7)sui7-L]; <Pa = 4>b +tan~ 1 [(—A/3/27) tan 7L]; and (f>B — A/3X/2. Using this 
representation, the coupling between the two waveguides can be regarded as a cascade of three processes: 1) 
phase retardation, 2) rotation, and 3) phase retardation. This becomes apparent if Eq. ^ is rewritten as 



exp(-j&,)T 3 T 2 Ti. 



with 



Ti 







cos (6/2) 
-j sin (6/2) 



-j sin (6/2) 

cos (6»/2) 



e -jT 2 





(3) 



(4) 



where Ti = 4>A — cf>B', ^2 = ~ < i>A~ ( t>B\ and Ti, T2, and T3 represent, in consecutive order, phase retardation, 
rotation, and phase retardation. The phase shift 4>b is a constant of no consequence. 

For perfect phase matching between the coupled modes, i.e., for A/3 = and an interaction coupling 
length L = qn/2K 1 where q is an odd positive integer, the coupling matrix T reduces to 



exp 



2 J 




-1 



-1 




(5) 



indicating that the modes are flipped. Applying this operation twice serves to double flip the vector, thereby 
reproducing the input, but with a phase shift twice that of qn/2. On the other hand, for 7L = pir, with p 



an integer, the matrix becomes 

cxp (jcj) A ) 



(-If 



exp(-j^) 



(6) 



Finally, for weak coupling (k « or k « A/3), we have 4>a ~ 0, whereupon T reduces to the identity matrix. 
Our interest is in three scenarios: 1) coupling between a pair of single-mode waveguides (SMWs); 2) 
coupling between a pair of two-mode waveguides (TMWs) ; and 3) coupling between a SMW and a TMW. 
The matrix described in Eq. @ is not adequate for describing the coupling in the latter two cases; in 
general, a 4 x 4 matrix is clearly required for describing the coupling between two TMWs. However, for the 
particular cases of interest here, the coupling between the two waveguides is such that only a single mode in 
each waveguide participates; this is because the phase-matching conditions between the interacting modes 
are either satisfied — or not satisfied. As an example for identical waveguides, similar modes couple whereas 
dissimilar modes fail to couple as a result of the large phase mismatch. The net result is that, for the cases 
at hand, the general matrix described in Eq. ^ reduces to submatrices of size 2x2, each characterizing the 
coupling between a pair of modes. 

5 Mode analyzer and modal Pauli spin operator a z 

A mode analyzer is a device that separates the even and odd components of an incoming state into two 
separate spatial paths. It is similar to the parity analyzer of one-photon parity space [40] , For the problem 
at hand, its operating principle is based on the selective coupling between adjacent waveguides of different 
widths. The even and odd modes of a TMW of width w\ are characterized by different propagation constants. 
An auxiliary SMW (with appropriate width W2, length L2, and separation distance b\ from the TMW) can 
be used to extract only the odd component [5] . The result is a mode analyzer that separates the components 
of the incoming state, delivering the the odd mode as an even distribution, as shown in Fig. Ufa). The 
end of the SMW is attached to an S'-bend waveguide, with initial and final widths W2, to obviate the 
possibility of further unwanted coupling to the TMW and to provide a well-separated output port for the 
extracted mode. If it is desired that the output be delivered as an odd distribution instead, another SMW 
to TMW coupling region (with the same parameters) may be arranged at the output end of the S'-bend, as 
illustrated in Fig. HJb) . This allows the propagating even mode in the SMW to couple to the odd mode of the 
second TMW, thereby delivering an odd distribution at the output. The appropriate coupler configuration 
is determined by the application at hand. It is important to note that the mode analyzer is a bidirectional 
device: it can be regarded as a mode combiner when operated in the reverse direction, as we will soon see. 

The Pauli spin (or spatial-parity) operator o~ z introduces a phase shift of 7r (imparts a negative sign) to 
the odd component of the photon state, leaving the even component unchanged; it thus acts as a half-wave 
retarder in mode space. It can be implemented by exploiting modal dispersion between the even and odd 
modes: a single TMW of length it/ \/3 e — /3 | , where /3 e and (3 are the propagation constants of the even and 
odd modes, respectively, results in the desired phase shift of n. For a weakly dispersive medium, however, a 
waveguide longer than practicable might be required. An alternative approach for implementing the Pauli 
spin operator o~ z involves cascading a mode analyzer and a mode combiner, as illustrated in Fig. HJc). As 
established in Eq. ([5]), perfect coupling between a pair of adjacent waveguides over an interaction length 
L = q-n /2k introduces a phase shift of qir/2, where q is an odd positive integer. A cascade of two such 
couplings thus results in a phase shift qn, with q odd, thereby implementing the Pauli spin operator a z . 
Proper design dictates that f3 e L e — (3 L , where L e and L are the distances traveled by the even and odd 
modes, respectively. Imperfections in the fabrication of the circuit may be compensated by making use of 
an electro-optic (EO) phase modulator, as sketched in Fig.[2jc). 

An example illustrating the operation of a mode analyzer, such as that shown in Fig.[2ja), is provided in 
Fig. [3] The behavior of the normalized propagation constants /3 of the even (to = 0) and odd (to = 1) modes 
before Ti indiffusion, as a function of the waveguide width u>, is presented in Fig. [3](a) for TM polarization 
at a wavelength of A = 0.812 /im. The horizontal dotted line crossing the two curves represents the phase- 
matching condition for an even and an odd mode in two waveguides of different widths. The simulation 
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Figure 2: (a) Sketch of a photonic circuit that serves as a mode analyzer (not to scale). It is implemented 
by bringing a single-mode waveguide (SMW) of width w-z and length Li into proximity with a two-mode 
waveguide (TMW) of width w\ . The two waveguides are separated by a distance b± . An S-bend waveguide 
of initial and final width W2, and bending length Lb, is attached to the end of the SMW. The center-to-center 
separation between the output of the S'-bend and the TMW is denoted S. All S-bends considered in this 
paper have dimensions Lb = 10 mm and S = 127/im (the standard spatial separation |26J. The odd mode is 
separated and delivered as an even distribution, (b) Sketch of a mode analyzer (not to scale) that separates 
the odd mode and delivers it as an odd distribution. It is more complex than the design presented in (a) 
because it incorporates a second TMW, again of width w\ , that is brought into proximity with a SMW of 
width u>2 and length L2 placed at the output of the S'-bend. These two waveguides are again separated by 
a distance b\ . (c) Sketch of a photonic circuit (not to scale) that changes the sign of the odd mode while 
leaving the even mode intact, thereby implementing the modal Pauli spin operator o~ z . An electro-optic 
phase modulator is used to compensate for any unintended differences in the phase delays encountered by 
the even and odd modes as they transit the circuit. 

presented in Fig. |3Kb) displays the evolution of the normalized amplitudes of the two interacting modes with 
distance. 



6 Mode rotator and modal Pauli spin operator a x 

The mode rotator is an operator that rotates the state by an angle 9 in mode space, just as a polarization 
rotator rotates the polarization state. It is also analogous to the parity rotator of one-photon spatial-parity 
space 40] . It achieves rotation by cascading a mode analyzer, a directional coupler, and a mode combiner; 
the three devices are regulated by separate EO phase modulators to which external voltages are applied. 
The mode analyzer splits the incoming one-photon state into its even and odd projections; the directional 
coupler mixes them; and the mode combiner recombines them into a single output. 

Implementation of the mode rotator is simplified by making use of the factorization property of the 
unitary matrix T that characterizes mode coupling in two adjacent waveguides (see Sec. 4). As shown in 
Eqs. © and Q, the coupling between two lossless waveguides can be regarded as a cascade of three stages: 
phase retardation, rotation, and phase retardation. If the phase-retardation components were eliminated, 
only pure rotation, characterized by the SO(2) operator, would remain. 

The phase-retardation components can indeed be compensated by making use of a pair of EO phase 
modulators to introduce phase shifts of T% and 1^, before and after the EO directional coupler, respectively. 
These simple U(l) transformations convert Ti and T3 in Eq. (TJ| into identity matrices, whereupon Eq. ((3J) 
becomes the SO(2) rotation operator. For a mode of wavelength A, and an EO phase modulator of length 
L and distance d between the electrodes, the voltage required to introduce a phase shift of T is V = 
XdT/Trrn 3 L , where the Pockels coefficient r assumes the values r 13 and r 33 , for n = n and n = n e , 
respectively [T71 Sec. 20. IB]. 

The standard EO directional coupler consists of two adjacent identical SMWs and makes use of an EO 
phase modulator to control the transfer of modal power between them p~7l Sec. 20. ID]. When no voltage 
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Figure 3: (a) Dependencies of the normalized propagation constants /3 of the fundamental (to = 0) and 
first-order (to = 1) modes on the widths w of the diffused channel Ti:LiNb03 waveguides. The input wave 
has wavelength A = 0.812 (im and TM polarization. The solid curves were obtained using the effective-index 
method described in [38], whereas the plus signs were computed using the software package RSoft. The 
dotted vertical lines represent the desired widths w\ and w 2 ■ (b) Simulated performance of a mode analyzer 
that takes the form displayed in Fig. [2ja). The blue curve represents the evolution with distance of the 
normalized amplitude of the odd mode in a TMW of width wi = 5.6 /zm, whereas the green curve shows 
the evolution of the even mode in a SMW of dimensions w 2 = 3.4 /im and L 2 = 6.2 mm. The separation 
between the TMW and the SMW is &i = 4//m and the iS-bend has dimensions L^ = 10 mm and S = 127 ^m. 
The dip in the curve for the SMW is associated with the tapered nature of the S'-bend. The results were 
obtained with the help of the software package RSoft. 
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Figure 4: Sketch of a photonic circuit that serves as a mode rotator (not to scale). It is implemented by 
sandwiching a directional coupler between a mode analyzer and a mode combiner. The coupling length 
of the directional coupler is ir/2n. To obtain a specified angle of rotation 9, voltages Vi, V 2 , and V3 are 
applied to the EO directional coupler, the input EO phase modulator, and the output EO phase modulator, 
respectively. 
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Figure 5: Operating voltages for the mode rotator vs. the angle of rotation 9. Voltages V\ (solid blue curve), 
V2 (dashed green curve), and V3 (dashed-dotted red curve) are applied to the EO directional coupler, the 
input EO phase modulator, and the output EO phase modulator, respectively. The input has wavelength 
A = 0.812 /im and TM polarization. The directional coupler comprises two identical SMWs separated by 
d = 5/xm; each SMW has width 2. 2 /mi and length 1.73 mm. The input and output EO phase modulators have 
electrode lengths of 5 mm and electrode separations of 5 /im. The curves represent theoretical calculations 
while the symbols represent simulated data obtained using the RSoft program. 

is applied to the EO modulator, the optical power is totally transferred from one waveguide to the other, 
provided that the interaction length L over which they interact is an odd integer multiple of the coupling 
length, it /2k [17j Sec. 8.5B]. The application of a voltage to the EO modulator introduces a phase mismatch 
between the two interacting modes that results in partial, rather than full, optical power transfer. In 
particular, if the voltage is chosen such that \A0L \ = \Z3tt (or V7n, \ZTTn, . . .), then no power is transferred 
between the two waveguides. The voltage required to introduce a phase mismatch of A/3 is approximately 
V = XdAf3/2Trrn 3 jTTJ Sec. 20. ID]. The waveguide beam combiner suggested by Buhl and Alferness [S3] 
operates on the same principle. 

However, because our modal state resides in a TMW, rather than in a SMW associated with the usual 
directional coupler, a mode analyzer with a configuration similar to that shown in Fig. Ufa) is used to direct 
the odd component to one arm of the EO directional coupler, and the even component to the other arm 
through an adiabatically tapered region, as shown in Fig. 01 A mirror-image tapered region and mode 
combiner follow the directional coupler to recombine the two components at the output of the device. 
Voltages Vi, V2, and V3 are applied to the EO directional coupler, the input EO phase modulator, and the 
output EO phase modulator, respectively. The voltages V2 and V3 can be modified as necessary to ensure 
that the overall phases acquired by the odd and even modes, both before and after the directional coupler, 
are identical when V\ = 0. 

An example showing the operating voltages V\ , V2 , and V3 required to obtain a specified angle of rotation 
9 is provided in Fig.[5j The directional-coupler voltage V\ has an initial value (for 9 — 0) that corresponds to 
a phase mismatch \A/3L | = y/3n; decreasing V\ results in increasing 9. When V\ = 0, the angle of rotation 
is 7r; the device then acts as the Pauli spin operator a x , which is a mode flipper (analogous to the parity 
flipper [40, 41 ). For V± = 0, there are an infinite number of solutions for the values of V% and V3, provided, 
however, that V2 = —V3. 

7 Controlled-NOT (CNOT) gate 

Deterministic quantum computation that involves several degrees-of-freedom of a single photon for encoding 
multiple qubits is not scalable inasmuch as it requires resources that grow exponentially |31j . Nevertheless, 
few-qubit quantum processing can be implemented by exploiting multiple-qubit encoding on single photons 
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Figure 6: Sketch of a Ti:LiNb03 photonic quantum circuit that behaves as a novel deterministic, two-qubit, 
single-photon, CNOT gate (not to scale). The control qubit is polarization and the target qubit is mode 
number. The circuit bears some similarity to the mode rotator shown in Fig. 2] both are implemented 
by sandwiching an EO directional coupler between a mode analyzer and a mode combiner. However, for 
the CNOT gate, the EO directional coupler comprises a pair of TMWs, whereas the mode rotator uses a 
conventional EO directional coupler utilizing a pair of SMWs. 



[54) . We propose a novel deterministic, two-qubit, single-photon, CNOT gate, implemented as a Ti:LiNb03 
photonic quantum circuit, in which the polarization and mode number of a single photon serve as the control 
and target qubits, respectively. 

The operation of this gate is implemented via a polarization-sensitive, two-mode, electro-optic directional 
coupler, comprising a pair of identical TMWs integrated with an electro-optic phase modulator, and sand- 
wiched between a mode analyzer and a mode combiner. It relies on the polarization sensitivity of the Pockels 
coefficients in LiNb03. A sketch of the circuit is provided in Fig. [6] The mode analyzer spatially separates 
the even and odd components of the state for a TM-polarized photon, sending the even component to one 
of the TMWs and the odd component to the other. At a certain value of the EO phase-modulator voltage, 
as explained below, the even and odd modes can exchange power. The modified even and odd components 
are then brought together by the mode combiner. 

To show that the device portrayed in Fig. [5] operates as a CNOT gate, we first demonstrate that the target 
qubit is indeed flipped by a TM-polarized control qubit, so that |1) = |TM). The polarization sensitivity of 
the Ti:LiNb03 TMWs resides in the values of their refractive indices n, which depend on the polarizations 
of the incident waves and the voltage applied to its EO phase modulator; and on their Pockels coefficients r, 
which depend on the polarization [T7| Example 20.2-1]. For a photon with TM polarization, the two-mode EO 
directional coupler offers two operating regions with markedly different properties. At low (or no) applied 
voltage, interaction and power transfer take place only between like-parity modes in the two waveguides 
because the propagation constants of the even and odd modes are different, so they are not phase-matched. 
However, at a particular higher value of the applied voltage, the behavior of the device changes in such a 
way that only the even mode in one waveguide, and the odd mode in the other, can interact and exchange 
power. This arises because the refractive indices of the two waveguides depend on the voltage applied to the 
device; they move in opposite directions as the voltage increases since the electric-field lines go downward in 
one waveguide and upward in the other. Figure [7] provides an example illustrating the dependencies of the 
propagation constants of the even and odd modes, in the two TMWs, as a function of the applied voltage. 

At a voltage indicated by the vertical dashed line in Fig. the even mode in one waveguide is phase- 
matched to the odd mode in the other. In a directional coupler with suitable parameters, a TM-polarized 
control bit will then result in a flip of the modal bit, whereupon cti\e) + a.2\o) — > a±\o} + a2|e). A TE- 
polarized control qubit, on the other hand, which sees n rather than n e , will leave the target qubit unchanged 
because of phase mismatch, so that |0) = |TE). Hence, the target qubit is flipped if and only if the control 
qubit is |1), and is left unchanged if the control qubit is |0), so that the device portrayed in Fig.[6]does indeed 
behave as a CNOT gate. In principle, it would also be possible to use a TE-polarized control qubit to flip 
the target bit; this option was not selected because it would require a higher value of EO phase- modulator 
voltage since the TE Pockels coefficient ri3 is smaller than the TM Pockels coefficient T33 [55] , 
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Figure 7: Dependencies of the normalized propagation constants f3 on the voltage applied to an EO TMW 
directional coupler comprising two waveguides [WG1 and WG2]. The propagation constants differ for the 
even and odd modes except at one particular voltage (vertical dashed line) where the even mode in one 
waveguide can be phase-matched to the odd mode in the other waveguide. The TMWs are identical, each of 
width 4 //m, and they are separated by 4 /mi. The input has wavelength A = 0.812 /mi and TM polarization. 
The symbols represent simulated data obtained using the RSoft program. 

A drawback of the photonic circuit illustrated in Fig. [6] is that it suffers from the effects of dispersion, 
which is deleterious to the operation of circuits used for many quantum information applications. Dispersion 
results from the dependence of the propagation constant /3 on frequency, mode number, and polarization. 
Polarization-mode dispersion generally outweighs the other contributions, especially in a birefringent material 
such as LiNbOs. 

Fortunately, however, it is possible to construct a photonic circuit in which the phase shifts introduced 
by dispersion can be equalized. A Ti:LiNb03 photonic quantum circuit that behaves as a novel dispersion- 
managed, deterministic, two-qubit, single-photon, CNOT gate is sketched in Fig. |SJ It makes use of three 
paths (upper, middle, and lower), in which the path- lengths of the three arms are carefully adjusted to 
allow for dispersion management. The third path provides the additional degree-of-frecdom that enables the 
optical path-lengths to be equalized. 

The design relies on the use of polarization-dependent mode analyzers at the input to the circuit. The TM- 
mode analyzer couples the odd-TM component of the state to the upper path, while the TE-mode analyzer 
couples the odd-TE component to the lower path. The even-TM and even-TE components continue along 
the middle path. Polarization-dependent mode combiners are used at the output of the circuit. 

If the control qubit is in a superposition state, the general quantum state at the input to the circuit, 
which resides in a 4D Hilbert space (2D for polarization and 2D for mode number), is expressed as 



|*i) = Qi|e,TM) + a 2 |o,TM) +a 3 |e,TE) + a 4 |o,TE) 
= |TM) <g> [ai|e> + a 2 |o) ] + |TE) ® [a 3 |e) + a 4 |o) ] 
= |e) <g) [ai|TM) + a 3 |TE) ] + \o) <g> [a 2 |TM) + a 4 |TE) ] , 



(7) 



where |e) and \o) are the basis states of the modal subspace; |TM) and |TE) are the basis states of the 
polarization subspace; the a's represent the basis weights; and (g) indicates the tensor product. Since the 
target (modal) qubit is flipped by a TM control qubit, the output state |<J/ ) becomes 



|* ) = ai|o,TM)+a 2 |e,TM)+a 3 |e,TE)+a 4 |o,TE) 
= |TM) ® [oi|o) + a 2 |e)] + |TE) <g> [a 3 |e) + a 4 |o) 



(8) 



where it is clear that the two terms in the input state, a i |e, TM) and a 2 \o, TM) , are converted to a ±\o, TM) 
and a 2 |e,TM), respectively, at the output, exemplifying the operation of this CNOT gate. Figure [3] displays 
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Figure 8: Sketch of a Ti:LiNb03 photonic quantum circuit that behaves as a novel dispersion-managed, 
deterministic, two-qubit, single-photon, CNOT gate (not to scale). The control qubit is polarization and the 
target qubit is mode number. The design is more complex than that shown in Fig.|6]because it accommodates 
dispersion management via path-length adjustments of the upper, middle, and lower paths. An EO TMW 
directional coupler is sandwiched between polarization-sensitive mode analyzers and polarization-sensitive 
mode combiners. The lower and upper waveguides of the two-mode directional coupler are denoted WG1 
and WG2, respectively. The paths taken by the components of the input state |^,) are shown, as is the 
output state |$ ). 

the paths taken by the components of the input state provided in Eq. ([7|); the output state set forth in 
Eq. (jHJ is also indicated. 

The output state in Eq. (jHJ is entangled in polarization and mode number; it is inseparable and cannot 
be written in factorizable form. A particular property of the CNOT gate is the induction of entanglement 
between factorized qubits: if the control qubit is in the superposition state -4=[ |TM) + |TE)], and the target 
qubit is in one of the computational basis states, then the output state of the CNOT gate is maximally 
entangled. An experimental test of the entanglement created between the polarization and modal degrees- 
of-freedom can be effected by using quantum-state tomography. The input to the CNOT gate can be readily 
generated from a product state, say |TM) ® |e), by rotation using a waveguide-based EO TE^TM mode 
converter [231155], in addition to a phase modulator, as described in Sec. 6. 

It remains to demonstrate the manner in which dispersion management can be achieved in the CNOT 
gate displayed in Fig. [5] The phase shift ip acquired by each component at the output is given by 



f e,TM = /3 e ,TM £\ + /3 ,TM ^2 + &' L D - {2qi + q 2 ) 7r/2 

f o,TM = <fi e,TM 

<p e ,TE =2/3 e ,TEh+/3"L D 

V o,TE = 2/3 0>T E h - 9371" + 2<t>A ! 



(9) 



where the /3's are the mode propagation constants; /?' is the propagation constant of either the TM-even 
mode in WG1 or the TM-odd mode in WG2; j3" is the propagation constant of the TE-even mode in WG1; 
<7i, <72, and q^ are odd positive integers that depend on the lengths of the TM-mode analyzer, directional 
coupler, and TE-mode analyzer, respectively; Lq is the length of the directional-coupler electrode; £\ is the 
path-length for the even modes before and after the directional coupler, £2 is the path-length for the odd- 
TM mode before and after the directional coupler; and 21 1 + Lp, 2£ 2 + Lr>, and 2^3 are the overall physical 
lengths of the middle, upper, and lower paths, respectively. The phase shift <pA arises from the coupling that 
affects the odd-TE component as it travels through the TM-mode analyzer. Phase shifts that accrue for the 
even modes as they pass through the mode analyzers and mode combiners are neglected because of large 
phase mismatches and weak coupling coefficients. By adjusting the lengths £\, £2, and £3, we can equalize 
the phase shifts encountered by each component of the state. Imperfections in the fabrication of the circuit 
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may be compensated by making use of EO phase modulators. 

A simulation that demonstrates the performance of the polarization-dependent mode analyzers and EO 
TMW directional coupler is presented in Fig- [HI The lengths £i, £2, and £3 are assumed to be adjusted such 
that they equalize the phase shifts encountered by each component of the state so that dispersion is not an 
issue. The spatial evolution of the normalized amplitudes of the odd and even modes inside the TM-mode 
analyzer, for TM- and TE-polarization, are displayed in Figs.^a) and[Ub), respectively. It is apparent that 
the TM-mode analyzer extracts only the TM-odd component, while the TE-odd component remains in the 
TMW waveguide until it couples to the lower path via the TE-mode analyzer [see Fig. M,o)]. Figures [SJd), 
(e), and (f) display the performance of the directional coupler for modal inputs that are TM-even, TM-odd, 
and TE-even, respectively. It is apparent in Fig.[§Jd) that the power in the even mode in WG1 is transferred 
to the odd mode in WG2 for TM polarization. Figure [^e) reveals complementary behavior: the power in 
the odd mode in WG2 is transferred to the even mode in WG1. Figure [9jf), on the other hand, shows that 
the TE-even mode travels through the directional coupler with essentially no interaction. Figures HJ^d), (e), 
and (f), taken together, along with the observation that the TE-odd mode preserves its modal profile during 
propagation, demonstrate a flip of the modal target qubit by the TM-polarized control qubit, and no flip by 
a TE-polarized control qubit, confirming that the photonic circuit in Fig. [5] behaves as a CNOT gate. 

The absence of a total power transfer from one waveguide to another in Figs. \§ld) and Ete) can be 
ascribed to sub-optimal simulation parameters. The conversion efficiency can be expected to improve upon: 
1) optimizing the length of the two-mode directional coupler; 2) minimizing bending losses by increasing 
the length of the S'-bend; 3) mitigating the residual phase mismatch by more careful adjustment of the 
voltage; and 4) improving numerical accuracy. Moreover, the deleterious effects of dc drift and temperature 
on the operating voltage and stability of the two-mode directional coupler can be minimized by biasing it 
via electronic feedback [56 ; a novel technique based on inverting the domain of one of its arms can also be 
used to reduce the required operating voltage [57] . Finally, it is worthy of note that decoherence associated 
with the use of a cascade of CNOT gates, such as might be encountered in carrying out certain quantum 
algorithms, may be mitigated by the use of either a qubit amplifier [32] or teleportation and error-correcting 
techniques [55] , 

8 Conclusion 

The modes of a single photon in a two-mode Ti:LiNb03 waveguide have been co-opted as basis states for 
representing the quantum state of the photon as a modal qubit. Various photonic quantum circuit designs 
have been presented for carrying out basic operations on modal qubits for quantum information processing 
applications. These include a mode analyzer, a mode rotator, and modal Pauli spin operators. We have also 
described the design of a deterministic, two-qubit, single-photon, CNOT gate, as well as a dispersion-managed 
version thereof, that rely on a single photon with both modal and polarization degrees-of- freedom in a joint 4D 
Hilbert space. The CNOT gate is a key element in certain sets of universal quantum logic gates. Simulations 
of the performance of all of these components, carried out with the help of the the commercial photonic and 
network design software package RSoft, provide support that they operate as intended. The design of these 
devices is based on selective and controllable power coupling among waveguides, the isomorphism between 
waveguide coupling and the SO(2) rotation matrix, and the tensor polarization properties of the Pockels 
coefficients in lithium niobate. The flexibility of Ti:LiNb03 as a material for the fabrication guided-wave 
structures should accommodate the development of increasingly complex quantum circuits and serve to foster 
new architectures. 
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Figure 9: Simulation demonstrating the performance of the polarization-dependent mode analyzers and 
the EO TMW directional coupler associated with the dispersion-managed, deterministic, two-qubit, single- 
photon, CNOT gate set forth in Fig. [8] The input wavelength is A = 0.812 /an. The TM-mode-analyzer and 
mode-combiner parameters are w\ = 5.6 /im, wi = 3.4 ^m, h\ = 4 ^m, and L2 = 6.2 mm; the TE-mode- 
analyzer and mode-combiner parameters are wi = 3/xm, b\ = 4/Ltm, and Li = 3.7 mm (see Fig. [2] for symbol 
definitions). The S'-bends have dimensions Lt = 10 mm and S = 127 /zm. The TMW directional-coupler 
has length L± = 2.2 mm, waveguide width w± = 5.6 /jm, electrode separation d = 4 /im, and an EO phase- 
modulator voltage V = 36 V applied to WG2, with WG1 at ground potential. All panels display the spatial 
evolution of the normalized amplitudes of the interacting modes, (a) The curves display strong coupling 
between the odd and even modes for TM-polarization inside the TM-mode analyzer. The input odd mode 
in the TMW is shown in blue and the even mode transferred to the SMW is shown in green [the same color 
conventions are used in panels (b) and (c)]. The even mode is ultimately coupled to another TMW at the 
output of the TM-mode analyzer and once again becomes odd. (b) The curves show negligible coupling 
between the odd and even modes for TE-polarization inside the TM-mode analyzer, (c) The curves display 
good coupling between the odd and even modes for TE-polarization inside the TE-mode analyzer. At the 
TE-mode combiner, the even mode in the SMW once again becomes an odd mode in the TMW. Panels (d), 
(e), and (f) display the performance of the directional coupler for modal inputs that are TM-even, TM-odd, 
and TE-even, respectively. For a given polarization, the blue and green curves represent the amplitudes of 
the even [denoted Even(l)] and odd [denoted Odd(l)] modes in WG1, respectively, while the the red and 
black curves are the amplitudes of the even [denoted Even(2)] and odd [denoted Odd(2)] modes in WG2, 
respectively. All simulated data in this figure were obtained using the RSoft program. 
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